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ABSTRACT 


This  work  reports  theoretical  investigations  of  the  surface  motion 
in  solids  caused  by  reflection  of  an  intense  acoustic  pulse  at  a  pressure 
release  surface  done  using  the  surface  velocity  reflection  coefficient. 
A  formula  was  derived  for  the  surface  velocity  reflection  coefficient 
which  relates  the  surface  particle  velocity  for  an  incident  compression 
wave  to  the  particle  velocities  of  the  incident  and  reflected  waves. 

To  obtain  the  surface  profiles,  two  different  assumptions  were  used. 
The  first  assumption  was  the  conventional  one.     In  this  assumption, 
surface  motion  was  neglected  during  the  period  that  the  pulse  strikes 
the  surface  and  it  is  called  the  fixed  surface  approximation.     The 
second  assumption  is  to  consider  the  motion  and  distortions  of  the 
surface  and  is  called  the  moving  surface  approximation. 

Four  different  materials  were  investigated  and  the  results  were 
compared  for  different  materials  and  in  the  two  different  approximations. 
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I.     INTRODUCTION 

Impulsive  loads  such  as  are  produced  by  explosions  and  impacts  may 
cause  the  failure  of  materials.     An  impulsive  load  is  distinguishable 
from  a  conventional  load  primarily  by  its  suddenness  of  application  and 
its  brevity  of  duration.     Because  of  these  loading  characteristics,  the 
reaction  of  a  solid  body  to  an  impulsive  load  must  be  understood  in 
terms  of  the  propagation  of  waves  in  the  body  generated  by  the  loading. 
If  the  loading  is  of  sufficiently  large  amplitude,  these  waves  produce 
extensive  fracturing  and  large  permanent  distortions  in  the  body  upon 
which  it  acts.     The  phenomena  associated  with  failure  under  impulsive 
loads  are  markedly  different  from  those  connected  with  failure  under 
conventional  loads.     The  reaction  of  materials  to  the  extremely  high, 
very  transitory  pressures  imposed  on  them  by  impulsive  loading  is  of 
great  fundamental  importance  [Ref.   1], 

The  basic  character  of  the  fracture  and  plastic  flow  which  are 
caused  by  impulsive  loading  is  beginning  to  be  understood  [Ref.   1]. 
A  common  fracture  process  results  by  the  generation  of  shear  waves 
whose  amplitude  exceeds  the  shear  strength  of  the  material  in  which 
the  v/ave  is  propagating.     Such  shear  waves  are  usually  generated  by 
the  reflection  of  a  strong  compressive  wave  from  boundaries.     Another 
common  fracture  process  results  when  the  pressure  release  wave  generated 
by  the  reflection  of  the  compressive  wave  from  a  boundary  has  too  great 
an  amplitude.     A  third  common  failure  process  involves  the  superposition 
of  shocks  reflected  from  various  boundaries.     Thus,  understanding  the 
way  in  which  intense  waves  damage  solids  requires  an  understanding  of 
reflection  of  waves  from  boundaries. 
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In  the  limited  theoretical  discussion  which  has  been  given  such 
processes  to  date  [Ref.  5],  the  reflection  process  is  discussed  using 
expression  obtained  for  small  amplitude  waves.     Experiments  show  that 
this  discussion  is  sufficient  for  qualitative,  but  not  always  quantita- 
tive, understanding  of  some,  but  not  all,  of  the  observed  phenomena. 

The  laws  governing  oblique  reflection  of  elastic  waves  at  a  free 
surface  are  well  known  [Refs.  2,  3].     There  is  a  particle  motion 
associated  with  the  reflection  of  these  elastic  waves  that  strike 
free  surfaces  obliquely.     The  object  of  this  thesis  is  to  study 
possible  departures  from  the  behavior  predicted  by  the  reflection 
formulae  which  result  because  of  the  finite  amplitudes  of  waves 
involved  by  investigating  the  surface  motion  in  solids  caused  by 
reflection  of  an  intense  acoustic  pulse  at  a  pressure  release 
surface. 


II.     SURFACE  MOTION  DUE  TO  REFLECTION  OF  ELASTIC  HAVES 

A.     ELASTICITY 

Mathematical   theories  of  elasticity  are  developed  in  the  same 
general  v/ay.     First,  the  notations  of  stress  and  strain  are  developed. 
Second,  a  stress-strain  relationship  between  these  quantities,  or 
their  derivatives,  is  assumed  which  idealizes  the  behavior  of  actual 
materials.     Finally,  using  this  relationship,  equations  of  motion  or 
equilibrium  are  set  up  which  enable  the  state  of  stress  or  strain  to 
be  calculated  when  a  body  is  subject  to  prescribed  forces. 

In  general,  the  stress  on  a  surface  element  in  a  solid  body 
does  not  act  normally  to  the  surface.     It  has  components  both  normal  to 
the  plane  and  tangential  to  it.     In  Fig.  2-1,  the  coordinate  system  has 
three  mutually  perpendicular  axes  0  ,  0  ,  0  .     The  stresses  acting  on 
three  planes  normal   to  these  axes  which  pass  through  a  point  P  are  con- 
sidered.    There  will  be  nine  components  of  stress.     These  are  denoted 
by  a     ,  a     ,  a     ,  a     ,  etc.     The  first  letter  in  the  suffix  indicates 
the  direction  of  the  stress  and  the  second  one  defines  the  plane 
in  which  stress  is  acting.     If  an  infinitesimal   rectangular  parallelepiped 
around  P,  with  its  faces  normal  to  the  axes,  is  considered,  by  taking 
moments,  it  may  be  seen  that  for  equilibrium  a      =  a     ,  a      =  a     ,  and 
a      =  a     .     Thus  there  are  only  six  independent  components  of  stress. 
These  completely  define  the  stress  at  the  point. 

The  strain  of  the  medium  is  the  relative  displacement  of  the  material 
particles  in  the  medium.     A  rigid  translation  or  rotation  of  a  medium 
is  not  considered  a  strain.     In  order  to  find  the  strain  at  the  point, 
its  change  in  position  relative  to  adjacent  points  must  be  analyzed. 
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Figure  2.1.  Stress  components  acting  on  an  infinitesimal  rectangular 

parallelepiped 
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Consider  a  point  very  close  to  P,  which  before  displacement  had 
coordinates  (x  +  dx),  (y  +  dy),  (z  +  dz).     Let  the  displacement  which 
it  has  undergone  have  components  (u  +  6u,  v  '+  6v,  w  +  &>),  where  u,  v 
and  (o  are  the  components  of  the  displacement  of  point  P.     If  6x,  6y 
and  6z  are  sufficiently  small. 


3X  3y     J        3Z 

(2.!)  5v  =  |£5x  +  H6y  +  ||62 

6w  =  3^  6x+ fy  6y+ 3z  6z 
The  relative  displacement  of  all  surrounding  points  P  may  be  found 

•  j-    ..  ,  £    ..  ....         3U      3U      3U      3V      3V      3V      3w      3w      3w 

if  the  values  of  the  nine  quantities  ^,  ^  ,p  ^  -^,  ^,  ■&'  ^,  ^ 

at  the  point  are  known.     These  nine  quantities  are  regrouped  as  follows: 


(2-2) 


=  iH. 

CXX   "     3X 

_    3V 

cyy  "  ay 

_   3io 
eZZ  *     3Z 

_  _3to    ,    3V 

eyz  "'  3y      3z 

,   3U    .    ou) 
£ZX   "    3Z         3X 

_    3V 
e        =  t —  + 

xy      3x 

3U 

ay 

2-      _   3w        3V 

x  "  3y  "  3z 

2"      _    3U        3^ 

wy  "  3z  "  3x 

2wz  -  ax  " 

3U 

sy 

The  first  three  quantities,  e     ,  e     ,  e__  correspond  to  the  fractional 

x  x      yy      z  z 

expansion.     The  second  three  e     ,  e        and  e      correspond  to  the  components 

of  shear  strain  in  the  planes  denoted  by  their  suffixes.     The  last  three 

a  ,  w    and  w    do  not  correspond  to  a  deformation  of  the  element  around  P, 
x      y  z 

but  are  the  components  of  its  rotation  as  a  rigid  body. 

Experimentally  it  is  found  that  for  most  solids,  the  observed  strains 
are  proportional  to  the  applied  load,  provided  that  the  load  does  not 
exceed  the  elastic  limit.     This  law  is  known  as  generalized  form  of 
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Hooke's  law,  and  can  be  stated  as:  Each  of  the  six  components  of  stress 
is  at  any  point  a  linear  function  of  the  six  components  of  strain.  This 
law  may  be  written  as 

a       -   Cue        +  Ct«e        +   Ci->e        +   C-i.e       +   C-,re       +   C-tz-e 
xx        llxx        12  yy        13  zz        14  yz        15  zx        16  xy 

ayy  ::  c21exx  +  c22Gyy  +  c23ezz  +  c24eyz  +  c25Ezx  +  c26£xy 


(2-3) 


zz         31   xx        32  yy         33  zz         34  yz         35  zx         36  xy 


yz        41  xx        42  yy         43  zz         44  yz        45  zx        46  xy 


zx        51  xx        52  yy         53  zz        54  yz        55  zx        56  xy 

axy  :    c61  xx  +  c62  yy  +  c63  zz  +  c64  yz  +  c65  zx  +  c66  xy 

The  coefficients  in  these  relations  are  known  as  the  elastic  constants 

of  the  material. 

In  an  isotropic  solid,   the  values  of  the  coefficients  are  independent 
of  the  set  of  rectangular  axes  chosen.     For  such  a  solid,  there  are  only  two 

independent  elastic  constants,  A  and  y  known  as  Lame's  constants.     All   the 

elastic  constants  can  be  expressed  in  terms  of  the  Lame's  constants. 

Thus,   the  Lame's  constants  completely  define  the  elastic  behavior  of  an 

isotropic  solid.     For  such  a  material,  equation  2-3  may  be  written 

°xx  =  XA  +  2lJExx      ayy  =  XA  +  2y£yy      azz  =  XA  +  2yczz 

(2-4) 

a       -   ye  a        -   ue  a        -   ue 

yz       M  yz  zx      M  zx  xy  xy 

Here  A  =  e      +  e       +  e     ,  and  is  called  dilatation.     It  represents  the 
xx         yy  z  z 

change  in  volume  of  a  unit  cube. 

For  convenience,   four  elastic  constants  are  normally  introduced. 
One  is  Young's  modulus  E,  which  is  given  by  a     /e       so  that: 


13 


(2-5)  z-.A^-LJA 

A  second  is  Poisson's  ratio  v  defined  as  -a     /e       so  that: 

yy    xx 

(2-6) 


7UT+ 

A  third  is  the  bulk  modulus  k  defined  as  the  ratio  between  the  applied 
pressure  and  the  fractional  change  in  volume  when  the  solid  is  subjected 
to  uniform  hydrostatic  compression  so  that: 

(2-7)  k  =  X  +  |H- 

Finally,  the  rigidity  of  shear  modulus,  which  is  identical  with  Lamp's 
constant  y,  is  defined  to  be  the  ratio  between  the  shear  stress  and  the 
shear  strain. 

B.  EQUATION  OF  MOTION  IN  AN  ELASTIC  MEDIUM 

In  order  to  obtain  the  equations  of  motion,  for  an  elastic  medium, 
the  variation  in  stress  across  a  small  parallelepiped  with  its  sides 
parallel  to  a  set  of  rectangular  axes  are  considered  (Fig.  2.2).  The 
components  of  stress  will  vary  across  the  faces;  to  obtain  the  force 
acting  on  each  face,  the  value  of  stress  at  the  center  of  each  face 
times  the  area  of  the  face  has  been  taken. 

Six  separate  forces  will  act  parallel  to  each  axis.  The  resultant 
force  acting  in  the  x-di  recti  on  is 

do 

(a       +  — **■  6x)6y&z  -  oY   <5y6z   + 
XX  9X  XX  J      „ 


+  KV  +   *lEL  <5y)$x6z  -   6      5X6Z  + 

Ay        3y  xy 
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"Z 


***£}>- 


Figure  2.2.     Stresses  acting  on  a  small  rectangular  parallelepiped 
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da 
+   (axz  +  -g~  6z)6x6y  -  axz6x6y 


3a  3a  3a 

v  3X  3y  3Z  y       J 


By  Newton's  second  law  of  motion,  neglecting  body  forces  such  as 
gravity,  this  will  be  equal  to 

(P«x6y6z)-^ 

where  p  is  the  density  of  the  element  and  u  is  the  displacement  in  the 
x-direction.     Thus  the  equations  of  motion  are 
2  3a  8a  3a 

f?    Q\  JO*    -  XX    +  W    +  XZ 

(2-8)  p-T  "  IT  +  IT  +  "TT 

3  t 

Similarly 

Jl  3a  3a  3a 

(2-9)  p^-4  =  -**  +  -^  +  -*£ 

,,2  3a  3a  3a 

/0   1M  3  u  _  zx   .   zy   ,        zz 

(2.10)  p-^  -  —  +  _JL  +  _ 

C.  ELASTIC  WAVES 

Equations  2-8,  2-9  and  2-10  will  hold  whatever  the  stress-strain 
behavior  of  the  medium.     By  substituting  the  equations  2-4  for  the 
stress  components  in  Equation  2-8,  one  obtains 

P7T  =  h  (XA  +  2yexx}  +  Jy  (yexy}  +  &  Kz> 

o  t 

By  using  the  equations  2-2,  one  finds 

(2-11)       p2-£=  (X   +  p)|f  +  yV2U 
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Similarly 

(2-12)  p^-4  =    (X  +   p)f£  +   yV2U 

(2-13)  pi-^=   (A  +  M)||+  yV2o) 

Equations  2-11   through  2-13  are  the  equations  of  motion  of  an 
isotropic,  elastic  solid  in  which  body  forces  are  absent.     They  predict 
the  propagation  of  two  types  of  waves  through  the  medium.     If  both 
sides  of  the  above  equations  are  differentiated  with  respect  to  x,  y  and 
z  respectively,  they  become 

2 
(2-14)  p^-4  =  (A  +  2p)v2A 

atr 

This  is  the  wave  equation  and  shows  that  the  dilatation  A  is  propagated 
through  the  medium  with  velocity 

(2-15)  c,  =  [(X  +  2u)/pf2 

On  the  other  hand,  if  A  is  eliminated,  it  can  be  shown  that 

(2-16)  p— **  =  yvV 

3t 

where  w    is  the  rotation  about  the  x-axis.     Similar  equations  may  be 

obtained  for  u>    and  Z  ,     Thus  the  rotation  is  propagated  with  velocity 


(2-17)  c2  =  (y/p)2 


This  shows  that  in  the  interior  of  any  elastic  solid  waves  may  be 
propagated  with  two  different  velocities.     These  two  types  of  waves 
should  be  termed  dilatation  waves  and  distortion  waves. 
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Any  plane  wave  propagated  through  on  isotropic  elastic  medium 
must  travel  one  or  the  other  of  the  above  velocities.     The  waves  which 
propagate  with  velocity  c-.  are  called  longitudinal  waves.     The  motion 
is  along  the  direction  of  propagation.     The  other  type  waves,  which 
propagate  with  velocity  c2,  are  called  transverse  waves  since  the 
motion  is  parallel   to  the  wave  front. 

The  equations  2-14  and  2-16  are  of  the  form 

(2-18)  A  m  C2V2^ 

at 

When  the  deformation  is  a  function  of  only  one  coordinate,  e.g.  x,  the 
equation  becomes: 


(2-19)  it  =  c2  it 

at  ax 


The  general  solution  of  this  is: 

ty  =  f(x  -  ct)  +  F(x  +  ct) 

where  f  and  F  are  arbitrary  functions  depending  on  initial   conditions, 
f  corresponding  to  a  plane  wave  travelling  along  the  positive  direction 
of  the  x-axis,  and  F  to  one  in  the  direction  opposite  to  this. 

D.     REFLECTION  AT  PRESSURE  RELEASE  SURFACE 

Normal  incidence.     Any  elastic  wave  will  be  reflected  when  it 
reaches  a  free  surface  of  the  material   in  which  it  is  travelling.     The 
simplest  case  occurs  when  the  wave  strikes  the  surface  normally.     In 
a  longitudinal  wave,  since  the  stress  normal   to  the  surface  at  the 
surface  must  be  zero,  the  reflected  pulse  must  be  opposite  in  sense 
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to  the  incident  pulse  (compression  reflected  as  tension  and  vice 
versa). 

Oblique  incidence.     When  an  elastic  wave  strikes  a  free  surface 
obliquely,  the  situation  is  much  more  complicated.     This  complication 
arises  from  trie  fact  that  the  energy  of  incident  wave  is  partitioned 
into  two  reflected  waves  instead  of  only  one.     When  either  a  longitudinal 
or  a  transverse  pulse  strikes  a  surface  obliquely,  both  a  longitudinal 
and  a  transverse  wave  are  reflected.     (Figure  2.3). 

Consider  first  a  longitudinal  wave  striking  a  surface  obliquely. 
There  will  be  a  reflected  longitudinal  portion  in  which  the  angle  of 
incidence  a  is  equal   to  the  angle  of  reflection.     There  will  also  be 
a  transverse  wave  which  is  generated  at  the  point  of  reflection. 
From  SnelVs   law,  the  relationship  between  the  angle  of  incidence  a 
and  the  angle  of  reflection  6  is 

(2"20)       Hrl =  § =  C2(1 " v)/(1 " Zv)fi 

From  equation  2-20,  it  is  seen  that  the  angle  of  reflection  6  is  the 
function  of  the  incidence  angle  and  the  Poisson's  ratio  of  the  material. 
Figure  2.4  shows  6  as  a  function  of  a  for  several  values  of  Poisson's 
ratio. 

This  generation  of  a  shear  wave  ca:         r  dualized  in  terms  of 
particle  velocities  (Fig.  2.5).     It  can  ajso-be  seen  that  a  shear 
wave  generated  by  the  oblique  reflection  of  a  longitudinal  wave  is 
polarized  in  the  plane  of  incident;  that  is,  the  particle  velocity 
associated  with  this  shear  wave  is  normal  to  the  direction  of  propaga- 
tion and  directed  in  only  one  direction.     The  case  of  a  compression 
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Figure  2.3    Oblique  incidence  of  a  plane  longitudinal  wave  at  a  free 

surface 
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wave  striking  the  free  surface  of  a  semi -infinite  body  will  be 
considered.     Let  a  be  the  instantaneous  value  of  normal  stress  of  the 
incident  wave,  a'   the  corresponding  normal  stress  of  the  reflected 
longitudinal  wave  and  t'   the  shear  stress  of  the  induced  transverse 
wave.     These  stresses  are  all  related  by  a  coefficient  of  reflection  R 
such  that  [Ref.  3] 

(2-21)  o'   =  Ra 

(2-22)  t'   =  [(R  +  1)   Cot2e]a 

and  the  coefficient  of  reflection  R  is  given  by   [Ref.  3]. 

2 

(2-2  3}  R  =  tan  B  tan  23  -  tan  a 

tan  3  tanfc23  +  tan  a 

R  is  a  function  of  a  and  3.     Since  R  is  a  function  of  3,  it  is  actually 
a  function  of  the  Poisson's  ratio  of  the  material.     Figure  2.6  shows 
R  as  a  function  of  a  for  several  values  of  Poisson's  ratios. 

A  very  similar  series  of  statements  can  be  made  for  the  oblique 
reflection  of  a  shear  wave.     By  referring  back  to  R  as  defined  by 
equation  2-23. 

(2-24)  t'  =  Rt 

(2-25)  a'  =  [(R  -  l)tan23]x 

where  t,  t'  and  a'  are  the  instantaneous  values  of  incident  shear 
stress,  generated  shear  stress  and  generated  longitudinal  stress, 
respectively,  3  now  representing  the  angle  of  incidence. 

When  an  incident  wave  strikes  a  free  surface  obliquely,  the 
surface  will  be  displaced  such  that  the  angle  of  emergence  o  satisfies 
[Ref.  3,  4]. 
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Figure  2.6.     Coefficient  of  Reflection  as  a  Function  of  Angle  of 
Incidence  for  Several  Values  of  Poisson's  Ratio 
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(2-26)  a  =  23 

The  angle  of  emergence  is  also  related  to  rectilinear  components 
of  the  particle  velocity  by 

(2-27)  a  =  tarf  1(vx/vy) 

where  v    is  the  sum  of  the  particle  velocity  components  perpendicular 

A 

to  the  surface,  v    is  the  sum  of  the  particle  velocity  vectors  parallel 
to  the  surface  as  seen  in  Fig.  2.5.     Again  it  can  be  seen  that  the 
angle  of  emergence  is  a  function  of  both  the  angle  of  incidence  and 
Poisson's  ratio. 

E.     SURFACE  VELOCITY  REFLECTION  COEFFICIENT 

The  surface  particle  velocity  for  an  incident  compression  wave  is 
related  to  the  particle  velocities  of  the  incident  and  reflected  waves. 
For  a  plane  elastic  longitudinal  wave,  the  particle  velocity  v*  at  any 
time  is  given  by  the  expression  [Ref.  4] 

(2-28)  vt  =  a/pc] 

where  a  represents  the  normal  stress  associated  with  the  wave.     For 
a  compressive  wave  the  particle  velocity  is  in  the  same  direction  as 
the  velocity  of  wave  propagation;  for  a  tension  wave  the  particle 
velocity  is  directed  oppositely  from  that  of  propagation  of  the  wave. 
A  corresponding  relationship  exists  for  particle  velocity  v.  within  a 
transverse  wave,  namely, 

(2-29)  vt  =  t/pc2 

where  t  is  the  shear  stress  associated  with  the  transverse  wave. 
From  Fig.  2.5,  the  surface  particle  velocity  is 
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(2-30)  S,  -  V  V  \ 

By  substituting  equations  2-28  and  2-29  into  the  equation  2-30 

where  n,  is  the  unit  vector  in  the  direction  parallel  to  the  incident 
longitudinal  wave,  n2  is  the  unit  vector  in  the  direction  parallel 
to  the  reflected  longitudinal  wave  and  ru  is  the  unit  vector  in  the 
direction  perpendicular  to  the  reflected  transverse  wave.     By  substitut- 
ing equations  2-21  and  2-22  into  the  equation  2-31 

(2-32)  vp  *  j2-  [n1  +  n2R  +  n3  J-  (R  +  l)Cot  2b] 


or 


(2-33)  vp  =  % 


where 


(2-34)  t  =  -~  [n,  +  n9R  +  n.  ■—  (R  +  l)Cot  2e] 


From  Fig.  2.5 


n,  =  l   Cos  a  +  j  Sin  a 
n«  ■  -l   Cos  a  +  j  Sin  a 
n3  =  i  Sin  6  +  j  Cos  3 

By  substituting  these  equations  into  equation  2-34,  it  can  be  shown 


that 
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(2-35)  t  =  4"(?[(1  '  R)  C°S  a  +  ^  °  +  R)  COt  23  S1n  3]  + 

pc1  L  2 

+  jftl  +  R)(Sin  a  +  J-  Cot  23  Cos  3)]| 

On  the  other  hand 

(2-36)  v    =  vp(i   Cos  26  +  j  Sin  23) 

By  using  equations  2-33,  2-35  and  2-36,  it  can  be  shown  that 

o       (1  +  R)    rCot  23  Cos  3  +  Sin  Bn  d 

(2-37)  R  =  i-rrr1-  ( SuTT3  }       a 


pc2 


V     =  Ra 
P 


where  R  is  called  the  surface  velocity  reflection  coefficient.     From 
this  formula,  it  can  be  seen  that  the  surface  velocity  reflection 
coefficient  is  a  function  of  the  angle  of  incidence,  the  Poisson's 
ratio  and  the  density  of  the  material. 
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III.     INVESTIGATION  OF  SURFACE  MOTION 

A.     NATURE  OF  INVESTIGATION 

The  theories  of  oblique  reflection  of  elastic  waves  at  a  free 
surface  were  outlined  in  the  previous  chapter.     Such  waves  cause  a 
particle  motion  on  the  surface.     The  motion  of  the  surface  will  cause 
the  reflection  characteristics  of  an  acoustic  pulse  to  change  as 
various  portions  of  the  pulse  encounter  to  surface.     The  theoretical 
investigations  of  the  surface  motion  in  solid  caused  by  reflection 
of  an  intense  acoustic  pulse  at  a  pressure  release  surface  were 
facilitated  by  using  the  surface  velocity  reflection  coefficient  R 
which  was  derived  in  equation  2-5. 

The  intent  was  to  illustrate  the  surface  motion  for  different 
materials,  for  different  loadings,  for  different  angles  of  incidence. 
For  this  purpose  the  IBM-360  computer  was  used.     The  surface  profiles 
were  obtained  by  using  two  different  approaches.     Tne  first  approach 
was  the  conventional  one  which  is  called  the  fixed  surface  approxima- 
tion.    This  assumes  that  surface  motion  can  be  neglected  during  the 
period  that  the  pulse  strikes  the  surface.     The  second  approach 
is  consider  the  motion  and  distortions  of  the  surface.     This  is 
called  the  moving  surface  approximation.     The  results  of  these  two 
approximations  were  then  compared. 

The  computer  programs  used  in  this  investigation  are  listed  in 
Appendix  A.     The  fixed  surface  approximation  computes  the  surface 
distortion  by  simply  integrating  the  surface  velocity,  obtained  by 
using  a  constant  surface  velocity  reflection  coefficient,  resulting 
from  a  pulse  arriving  at  the  x  =  0  plane.     The  moving  surface  approxima- 
tion makes  two  major  corrections.     First,  the  surface  velocity  is 
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computed  at  each  time  from  the  displaced  position  of  the  surface, 
the  displacement  resulting  from  the  previous  portions  of  the  pulse. 
Secondly,  account  was  taken  of  changes  in  the  angle  of  incidence 
caused  by  distortions  of  the  surface  from  a  plane.     This  was  done 
by  constructing  the  normal   to  an  equation  describing  the  surface 
obtained  by  passing  a  quadratic  curve  through  the  point  in  question 
and  its  neighboring  points.     Because  of  the  time  and  position 
dependent  changes  in  the  pulse  shape,  numerical  integration  was 
used  to  convert  the  surface  velocity  into  surface  displacement. 

Throughout  tin's  investigation  a  saw-toothed  shape  stress  pulse 
with  1  mm.  pulse  width  has  been  used  to  illustrate  the  surface  motion. 
The  amplitude  of  the  stress  was  changed  from  a  value  equal   to  the  normal 
fracture  stress  to  ten  times  this  quantity.     The  normal  fracture  stress 
is  the  maximum  tensile  stress  which  the  material  will   tolerate.     A  more 
exact  definition  would  be  the  minimum  dynamic  tensional  stress  required 
to  rupture  the  material.     The  angle  of  indidence  was  varied  from  5  to 
85  degrees,  using  5-degree  increments.     Four  materials,  aluminum,  steel, 
plexiglass,  and  lead,  were  investigated.     The  characteristics  of  these 
materials  are  given  in  Table  I. 

B.     RESULTS 

Comparison  of  Approximations:     As  an  illustration  of  the  differences 
between  the  two  approximations,  the  surface  motion  produced  in  aluminum 
and  lead  for  a  sawtooth  compressive  pulse  of  amplitude  equal  to  the 
fracture  stress  incident  at  30°  computed  in  the  two  approximations  is 
illustrated  in  Figs.   3.1  and  3.2.     The  maximum  displacement  produced 
by  the  pulse  is  the  same  in  the  two  approximations.     This  result  is 
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Figure  3.1.     Surface  motions  for  aluminum  under  two  different  assumptions. 
Amplitude  of  stress  is  equal   to  normal   fracture  stress,  a  =  30' 
time  =  0.251  x  10-6  sec. 
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Figure  3.2.  Surface  motions  for  lead  under  two  different  assumptions. 
Amplitude  of  stress  is  equal  to  normal  fracture  stress,  a 
time  =  0.741  x  10"6 


=   30°, 


32 


expected  since  the  net  displacement  can  be  computed  from  the  impulse 
contained  within  the  incident  pulse  and  the  elastic  moduli  of  the 
medium  within  which  it  propagates.     The  moving  surface  approximation 
predicts  a  more  rapid  displacement  of  the  surface  than  the  fixed 
surface  approximation.     This  also  is  quite  understandable. 

The  Angle  <j>:     The  change  of  the  angle  of  incidence  due  to  surface 
distortion  is  called  $,     During  the  motion  of  surface  <}>  at  any  given 
spot  changes  with  time.     To  illustrate  the  change  in  <t>  by  the  time  at 
a  given  position,  the  curve  is  drawn  for  aluminum  and  given  in  Fig.  3.3. 
It  shows  that  the  angle  $  starts  to  increase  after  the  pulse  strikes 
to  the  surface,  and  reaches  its  maximum  value,  then  decreases  again. 
But  the  decrease  is  relatively  slower  than  the  increase.     This  is  a 
result  of  the  assumed  sawtooth  pulse.     The  angle  <j>  is  also  effected 
by  the  angle  of  incidence  and  by  the  amplitude  of  stress.     The 
effect  of  the  angle  of  incidence  is  illustrated  in  Fig.  3.4. 
Figure  3.4  shows  the  change  in  angle  <j>  by  the  angle  of  incidence  for 
all  materials  investigated.     The  angle  $  changes  linearly  with  stress 
amplitude  because  of  using  the  Hooke's  Law  and  the  relatively  small 
displacements.     The  change  in  the  angle  <}>  with  stress  amplitude  is 
given  in  Fig.   3.5.     In  this  figure,  aluminum  is  chosen  as  an  example. 

In  the  moving  surface  assumption,  the  distortion  of  the  surface 
from  a  plane  effects  the  pulse  shape  and  causes  distortion  on  the  pulse. 
Assuming  plane  wave  propagation,  there  is  no  distortion  at  the  front 
and  tail  end  of  the  reflected  pulse  since  it  is  reflected  from  a  flat 
surface.     But  in  between  these  two  ends,  the  different  portions  of  the 
pulse  cause  different  distortions  on  the  surface  which  result  in 
different  angle  of  correction  <j>.     Because  of  different  $,  the  reflected 


33 


- 


& 


O) 

en 

i 

o 

'x 

-e- 


-1 


-2 


-3 


-4 


time  (xlO    sec) 

*".  *      1 

322 

r 

j^ij 

Figure  3.3.     Change  in  angle  <f> with  time  for  aluminum.     Stress  amplitude 
equals  to  normal  fracture  stress. 
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angle  is  not  the  same  as  the  instant  angle  of  incidence  for  these 
portions  of  pulse.     Thus  the  wave  fronts  for  these  portions  of  pulse 
are  not  in  the  same  direction  relative  to  the  front  and  tail  end  of 
the  pulse.     Hence  the  reflected  pulse  under  the  moving  surface 
assumption  is  not  the  same  as  the  incident  pulse.     From  the  investiga- 
tion, it  was  seen  that  the  pulse  became  far  apart  after  the  reflection 
from  the  surface.     And  this  causes  the  distortion  in  the  stress  field 
which  is  not  investigated  in  here. 

The  moving  surface  assumption  also  effects  the  amplitude  of  the 
reflected  pulse.     The  amplitude  of  the  reflected  stress  pulse  equals 
to  the  product  of  the  coefficient  of  reflection  and  the  amplitude  of 
the  incident  stress  pulse.     But  from  Fig.  2.6,  it  can  be  seen  that 
the  coefficient  of  reflection  is  the  function  of  the  angle  of  incidence. 
Since  the  assumption  causes  a  change  in  the  angle  of  incidence,  the 
amplitude  of  the  reflected  pulse  is  not  the  same  as  the  amplitude  of 
the  reflected  stress  pulse  under  the  fixed  surface  assumption  condition. 
It  is  either  weaker  or  stronger  than  it.     It  depends  on  the  angle 
of  incidence  and  the  material  used.     To  illustrate  the  effect  of  <j> 
on  the  amplitude  of  stress  pulse,  two  numerical  examples  were  done 
for  aluminum.     The  coefficients  of  reflection  were  calculated  by 
using  the  equation  2-23  for  two  different  angles  of  incidence,  30 
and  40  degrees.     The  amplitude  of  stress  was  the  normal  fracture 
stress  and  the  maximum  angles  of  correction  were  used  in  both  cases. 
For  30  degrees  angle  of  indidence,  by  using  fixed  surface  approximation, 
the  coefficient  of  reflection  equals  -0.778.     But  in  the  same  condition 
by  using  the  moving  surface  approximation,  it  equals  -0.790.     For 
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45  degrees  angle  of  incidence,  the  first  approximation  results  in  the 
coefficient  of  reflection  equals  to  -0.579,  the  second  approximation 
results  in  the  coefficient  of  reflection  equals  to  -0.590. 

Surface  Displacement:     Surface  displacement  is  effected  by  the 
stress  amplitude  and  by  the  angle  of  incidence.     The  surface  displace- 
ment increases  linearly  with  stress  amplitude.     However,  it  decreases 
by  increasing  the  angle  of  incidence.     The  surface  motions  for  two 
different  materials,   aluminum  and  lead,  are  illustrated  in  Figs.  3.6 
and  3.7. 
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Surface  motions  for  aluminum.  Stress  amplitude  equals  to 
normal  fracture  stress  and  a  =  30. 
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APPENDIX  A 


THE  COMPUTER  PROGRAMS 

FOR 

SIMULATION  OF  SURFACE  MOTIONS 
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